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Abstract. Supersymmetry (SUSY) is an attractive extension of the Standard Model (SM) of particle 
physics which solves the SM hierarchy problem. Motivated by the theoretical /x-term problem of the 
Minimal Supersymmetric Model (MSSM), the Next-to MSSM (NMSSM) can also account for experimental 
deviations from the SM like the anomalous muon magnetic moment and the dark matter relic density. 
Natural SUSY, motivated by naturalness considerations, exhibits small fine tuning and a characteristic 
phenomenology with light higgsinos, stops and gluinos. We describe a scan in NMSSM parameter space 
motivated by Natural SUSY and guided by the phenomenology of an NMSSM with a slightly broken 
Peccei-Quinn symmetry and a lightly coupled singlet. We identify a scenario which survives experimental 
constraints with a light singlet Higgs and a singlino lightest SUSY particle. We then discuss how the 
scenario is not presently excluded by searches at the Large Hadron Collider (LHC) and which channels are 
promising for discovery at the LHC and International Linear Collider. 

PACS. 11.30.Pb Supersymmetry - 14.80.Da Supersymmetric Higgs bosons - 12.60.Jv Supersymmetric 
Models 


1 Introduction 

With the discovery of the 125 GeV Higgs boson h \25 by 
ATLAS |lj and CMS [2] at the Large Hadron Collider 
(LHC), particle physics enters a new era. In the Standard 
Model (SM) of particle physics, the properties of the Higgs 
boson are determined by theory once the mass is known 
|3]. At present, their measurements are consistent with the 

But the SM is not complete. Experimentally, it does 
not account for Dark Matter (DM), the anomalous muon 
magnetic moment or the strong CP problem, among other 
things. Theoretically, it suffers from the hierarchy prob¬ 
lem. Supersymmetry (SUSY) solves the hierarchy prob¬ 
lem by introducing a fermionic partner for each SM boson 
and a bosonic partner for each SM fermion [HD . SUSY 
with conserved R parity provides a natural candidate for 
DM, the Lightest Supersymmetric Partner (LSP), and can 
account for the anomalous muon magnetic moment by in¬ 
troducing new particles in loops. 

The principle of naturalness in physics maintains that 
an effective physical theory approximately valid below some 
characteristic scale should not be very sensitive to the cor¬ 
rect theory above that scale m- Applied to electroweak 
symmetry breaking in SUSY, this implies that the success 
of the effective SM Higgs theory disallows SUSY too far 
above the electroweak scale m- In particular, the char¬ 
acteristic mass spectrum of Natural SUSY includes light 
superpartners of the Higgs bosons, top quark and gluon 
near the electroweak scale. 


The Minimal SUSY Model (MSSM) contains only the 
SM particles and their superpartners, together with an en¬ 
larged Higgs sector: one neutral pseudoscalar, two neutral 
scalars and two charged scalars which arise from the two 
Higgs doublets H u and Hd necessary for the Higgs mecha¬ 
nism in SUSY [3]. But the MSSM suffers from the so-called 
/i-term problem, which prevents the term [iH u Hd in the 
MSSM superpotential from reaching the electroweak scale 
without fine tuning mmm- 

The Next-to MSSM (NMSSM) solves the /x-term prob¬ 
lem by introducing a singlet S and replacing [iH u Hd with 
A SH u Hd- The Z 3 invariant NMSSM superpotential is [T31 

M 

W = XSH u H d + | S 3 (1) 

where A and k are free parameters. An effective /x-term is 
generated as the vacuum expectation value of S , /x e // = 
A (S'), reaching a natural scale without fine tuning |141I13| . 

In addition to the Higgs content of the MSSM, the 
NMSSM contains an additional pseudoscalar and an ad¬ 
ditional scalar so that the NMSSM Higgs sector consists 
of two neutral pseudoscalars (ai, 02 ), three neutral scalars 
(hi,/i 2 ,h 3 ) and two charged scalars ( H + ,H ~) [TMTfi] . 
The NMSSM Higgs sector is fully determined at tree level 
by A and k , A\ and A K (soft trilinear couplings), /x e // and 
tan/3 (ratio of H u ,Hd vacuum expectation values) [15]. 

One notable version of the NMSSM is the Peccei-Quinn 
(PQ) symmetric NMSSM, characterized by k = 0 [TT11T51 






2 


C.T. Potter: Natural NMSSM with a Light Singlet Higgs and Singlino LSP 


mm■ The PQ symmetric NMSSM explains why there is 
so little CP violation in the strong sector by exhibiting an 
axion, the massless pseudoscalar ai. In the NMSSM with 
a slightly broken PQ symmetry, with small k and A K , the 
ai acquires a small mass proportional to kA k but can still 
solve the strong CP problem P~51fT51[T5LI20L[^T] . 

Scenarios with a light NMSSM pseudoscalar Higgs, 
motivated variously by the strong CP problem, natural¬ 
ness, the anomalous muon magnetic moment, the mass 
spectrum, and the similarity of the baryon density to the 
dark matter density, have been discussed in the literature 
n^^ll^llMIf^llTfi 1127 1I2%1 . In this study we assume a light 
NMSSM pseudoscalar ai with 2 m T < m ai < 2mg. Moti¬ 
vated by the LEP Zbb feature near m bb ~ 60 GeV [25], we 
identify this as an hi candidate. We further identify the 
ft. 125 as the second lightest neutral scalar /12 of the NMSSM 
and note that the fti 25 signal strength measurements at 
the LHC [6]|5j place the heavier NMSSM 02 , h 3 , H + in the 
effective MSSM decoupling limit. 


2 Effective MSSM (A, k ps 0) 

We now consider the phenomenology of the NMSSM with 
a slightly broken PQ symmetry in which the singlet S is 
completely decoupled from the doublets H u and Hd (A = 
0). We then consider how the phenomenology is altered 
when the singlet is allowed a weak coupling to the doublets 
(A ps 0). The case A, k ps 0 is known as the effective MSSM 

m- 

For A = 0, there is no mixing of the singlet with the 
doublets. The generic couplings in the NMSSM have been 
detailed in mm- We adopt the notation of the latter, 
denoting S?- (P?) as the jth component of mass eigenstate 
hi (a^, where j = 1 , 2 ,3 corresponds to the u doublet, the 
d doublet, and singlet respectively. For purely singlet fti 
and ai, £>13 = Pi 3 = 1 and all other Sij, Pij vanish, so 
the ai cannot decay to SM particles since their coupling 
is proportional to Pn = 0 or Pi 2 = 0 , and similarly for 
the fti. The ai is stable and the only allowed fti decay 
for mfcj ~ 60 GeV is fti -» oidi. The singlet sector is 
decoupled from the SM sector. 

Furthermore, for the case A = 0, the singlet sector is 
decoupled from the MSSM sector. One neutralino is pure 
singlino whose mass, at tree level, is related to the ai 
mass by m x = —2m l ax l?>A Kj [T71H5] , For m ai ~ 10 GeV 
and \A k \ of 0(1) GeV, consistent with a slightly broken 
PQ symmetry, this yields m x ps 60 GeV. In this study we 
identify the singlino as the LSP xi- Denoting Nfj as the 
jth component of Xi, where j = 1 , 2 ,3,4,5 corresponds 
to bino, wino, u higgsino, d higgsino, and singlino respec¬ 
tively. Neutralinos heavier than the singlino LSP have zero 
singlino component, JV 15 = 1 and all other /V )5 vanish. No 
heavier neutralino can decay to the singlino since the cou¬ 
pling is proportional to N i5 = 0 for i > 1. The NLSP X 2 
is stable for conserved R parity. 

However, when the singlet is allowed a weak coupling 
to the doublets (A ps 0), small mixing between the sin¬ 
glet sector and the SM and MSSM sectors is possible. 


In this case the a\ may couple to SM pairs so the a\ is 
no longer stable. For m ai ps 10 GeV, ai —> r + r” dom¬ 
inates, with decays to gluon and light quark pairs sub¬ 
dominant. For mftj ps 60 GeV, fti — > a 3 ai remains domi¬ 
nant with decays to SM pairs, notably fti —> bb , subdom¬ 
inant. Furthermore, the xi can be produced from heav¬ 
ier neutralino decays since singlino mixing with doubli- 
nos is allowed. Then the NLSP X 2 and NNLSP x '3 are 
no longer stable against decay to Xi- Above the threshold 
m x2 = m Xl +m ai ps 70 GeV, X 2 —> Xi a i dominates while 
below it X 2 —> XiZ* dominates. In the latter case, the de¬ 
cay may occur outside the LHC detector effective tracking 
volume for A of 0(1O -2 ) or less [31] . Above the threshold 
m x 3 = m Xl + TOftj ps 120 GeV, X 3 —t Xi^i dominates 
while below it X 3 —> Xi a i an< 4 X 3 ~t Xi^* dominate. 

Further information about A and n can be extracted 
from the ft] .2 sum rule m 

+ ml 2 ps to| + ^kv s (4kv s + V2A K ) (2) 

where v s = s/2[x e f //A. For nq,, = 60 GeV, m,h 2 = 125 GeV 
and n e ff = 300 GeV, the sum rule yields |k/A| ps 0.176. 

To summarize, we assume an effective MSSM with 
mostly singlino LSP xi an d m Xl ps 60 GeV. The a\ and 
fti are mostly singlet with dominant decays to SM r pairs 
and a\ pairs, respectively. The a±, hi and xi can be pro¬ 
duced in neutralino decays. For m X2 ps 70 GeV or be¬ 
low and A < 0(1O~ 2 ), the X 2 decays outside of the ef¬ 
fective tracking volume. For m X3 ps 120 GeV or above, 
X'3 —> X 1 A 1 i s dominant. Finally, the fti^ mass sum rule 
yields |«/A| ps 0.176 for fi e ff = 300 GeV. These consid¬ 
erations, together with naturalness, inform the parameter 
ranges in the scan described in the next section. 


3 Parameter Scan 

The parameter scan is performed with NMSSMTools4.4.0 
CT21[T{Tnimi35II36IfTfl . probing 10 8 random points. We trade 
the soft trilinear parameters A\,A K ,A t for mp,mA,X t , 
defined by mm 


sin2/] (^ + '».) 

(3) 

3 1 
-~j=K v sA k 

(4) 

A t - n e ff/ tan/3 

(5) 


Here (mp) is the diagonal component of the CP odd 
doublet (singlet) mass matrix and X t is the stop mixing 
parameter. 

The parameters scanned are A, k, rriA, mp, n e ff, tan/3, 
M 2 , X t and toq 3 . We fix the gaugino masses Mi and M 3 
with the unification constraints Mi = 5 M 2 and M 3 = 
3M 2 . We further assume 77 iq 3 = mp 3 - All other squark 
and soft trilinear parameters are fixed to 1500 GeV, and 
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Parameter 

Range/Value 

h60 

A 

(0,0.1] 

0.03505 

K , 

[-0.01,0.01] 

0.006088 

m A 

[500,1500] GeV 

1068. GeV 

mp 

[9.9,10.5] GeV 

10.25 GeV 

Me// 

[100,300] GeV 

166.7 GeV 

tan/3 

[1,30] 

15.49 

Mi 

iM 2 

80.73 GeV 

m 2 

[100,300] GeV 

161.5 GeV 

M 3 

3 M 2 

484.4 GeV 

x t 

[0.8X? ax , 1.8Xr x ] 

1378. GeV 

m Q 3 L 

[350,550] GeV 

546.9 GeV 

rn TJS R 

m Q3 

546.9 GeV 


Table 1 . NMSSM parameters and their scan ranges. Addi¬ 
tionally, k is constrained to satisfy 0.125A < |k| < 0.225A. The 
point heo (k = 0.006088 and A K = —1.087 GeV) is taken from 
points surviving the scan and is described in Section [o] 


the slepton mass parameters are fixed to 200 GeV. See 
Table Q] for scanned parameter ranges. 

Motivated by the PQ symmetric NMSSM, we scan 
small k and A K or equivalently, from Equation [I] small k 
and small mp. The lower range bound of mp (9.9 GeV) is 
informed by the anomalous muon magnetic moment study 
[25] , while the upper bound (10.5 GeV) is informed by the 
r]b mass spectrum study m- Then k is scanned in the 
range — 0.01 < k < 0.01 and is also required to satisfy 
0.125A < |re| < 0.225A since this requires ~ 60 GeV 
within several GeV. We scan moderately small A in the 
range 0 < A < 0.1. Since the /ii 25 signal strength con¬ 
straints are applied in the scan, tua is allowed to go into 
the effective MSSM decoupling limit ipa rriz to accom¬ 
modate the SM-like couplings of the / 1125 . 

The neutralino and chargino masses are largely deter¬ 
mined by fi e f /, M\ and M 2 which, from naturalness con¬ 
siderations are bounded above in the scan by 300 GeV 
[ 12 ]. At tree level, the stop masses are m| ^ = mg 3 + 

m 2 ± mtX t for toq 3 = mu 3 [15l |. Naturalness informs the 
rriQ 3 range since light stops are compatible with small fine 
tuning. 

The tree level Higgs mass in the MSSM is bounded 
by to 2 < m? z cos 2 2/3, requiring a large loop correction for 
the /ii 25 - In the NMSSM the upper bound on m 2 has an 
additional 0(X 2 v 2 ) term. The stop mixing parameter X f 
partly determines the one loop correction m- 


Sm 2 h 


3 G F 
•\/2 7T 2 


m t 





where the parameter toj is defined by m~ = \(rn 2 ^ +m | ). 
The correction is strongly dependent on the top mass mt- 
In the scan m* = 172.5 GeV. 

To allow the large correction required by the / 1125 , but 
with small mf required by Natural SUSY, the stop mix¬ 
ing X t is allowed to contribute up to its maximal possible 
correction at X™ ax = In the scan NMSSMTools4 


calculates the Higgs mass spectrum at one-loop level in¬ 
cluding external momentum for self-energies and two-loop 
level excluding external momentum | i35II39j . 


4 Surviving Points 

The suite of constraints imposed by NMSSMTools4 while 
scanning includes experimental results from a wide variety 
of sources, including: 

— Anomalous muon magnetic moment Aa M measured by 
BNL E821 gO] 

— DM relic density measured by Planck [4T], di¬ 

rect DM exclusion by LUX m 

— B Physics, b —> sj, B —> X s /r + /i , B + —> t + p, B s —> 

T(ls) -> 07 , ??b(ls) 

— Higgs. LHC hi 25 , LEP e + e _ — Zh, Tevatron/LHC 
t bH + , NMSSM searches 

— SUSY. Tevatron/LHC y + , <Z, < 7 , e, fi, f mass constraints, 
t —)• blv, x° c an d b —)• x°b 

Loose constraints imposed during the scan require rrq , 2 a: 
125 GeV within 3 GeV and impose an upper bound on 
each hi 25 signal strength y 2 , calculated as in [43] . Of the 
10 8 points scanned, 42 survive the constraints imposed 
during the scan. 

Constraints are tightened after the scan. The low mass 
Higgs sector must satisfy: 

9.9 < m ai < 10.5 GeV 
50 < m hl < 70 GeV 
122 < m h2 < 128 GeV 

Finally, the sum of /1125 signal strength y 2 are required to 
satisfy )T/ y? < 13. Of the 42 points surviving the scan 
constraints, 15 points survive these final constraints. 

In order to demonstrate the naturalness of the sur¬ 
viving points, we examine the fine tuning metric F max = 

maXaeA ( ) emulated by NMSSMTools4. This met¬ 

ric yields the largest fine tuning over fundamental pa¬ 
rameter set A. Surviving points have small fine tuning, 

5 < F max < 10, light stops 300 < < 400 GeV and 

light gluinos 500 < rrig < 650 GeV. While agreement is 
not universal on which F max values characterize low fine 
tuning m , studies have considered F max of order 0(1O 2 ) 
to be typical for the NMSSM [44] and (^(lO 1 ) to be low 
fine tuning |24ll23ll22ll45] . A recent study seeking to es¬ 
tablish naturalness as objective, model-independent and 
predictive concludes that a SUSY model with F max < 30 
is natural, while one with F max < 10 is stringently natural 

US- 


5 Benchmark h 6 o 

It has been noted that points surviving the scan repre¬ 
sent a Natural NMSSM with slightly broken PQ sym¬ 
metry. They also exhibit a light pseudoscalar Higgs with 
m Ql « 10 GeV, a light scalar Higgs with m ^ « 60 GeV, a 
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NMSSMTools 4.4.0 
PySLHA 3.0.2 



Fig. 1. Masses and decays for the lighter heo benchmark spec¬ 
trum. Only decays with branching ratios larger than 5% are 
shown. Note that the /12 is mostly decoupled. 


Component 

Ol 

hi 

Xi 

X2 

X3 

Pu/Sn/Nn 

-0.002 

-0.006 

0.151 

0.882 

0.401 

Pn/Sn/Na 

0.000 

-0.130 

-0.054 

-0.153 

0.679 

P 13 /Si 3 /N i5 

1.000 

0.992 

0.981 

-0.189 

0.045 


Table 2. Doublet and singlet components of the ai, hi and 
gaugino and singlino components of the XiiX2,X3 > n ^ 60 - 


singlino LSP DM candidate with m Xl ss 60 GeV annihilat¬ 
ing via X 1 X 1 -;► bb, and a light stop with ~ 350 GeV. 

The benchmark point heo satishes the threshold crite¬ 
rion m x3 > + m Xl with the largest branching ratio 

for X‘-i —>• Xi h\ of all surviving points in the scan. The low¬ 
est branching ratio for this decay in the surviving points 
which reach threshold is 65%, while the highest is 80%. 
This ensures production of ai from hi —> oiai in stop 
pair events with t\ —> xtb —► X3 Wb or t\ —»• \ 3 t. See 
the last column of Table Q] for the numerical values of 
the parameters which define heo- See Figure [1] generated 
with PySLHA [27| using the SUSY Les Houches Accord 
(SLHA) pf8lRH7] file produced by NMSSMTools4, for the 
mass spectrum and decays in heo- 

For components Pij, Sij, Nij of the ffli,/ii,XijX 2 and 
X 3 in heo see Table EJ The LSP xi is dominantly singlino, 
while the a\ and h\ are dominantly singlet. The Xi mixing 
with doublinos and gauginos is small, as is the mixing of 
the ai and hi with the doublets. The phenomenology of 
a singlino LSP at the LHC has been considered in EDE23 
The phenomenology of light stops 
in the NMSSM, and how they avoid exclusion at the LHC, 
has been recently considered in [55] , 

For the numerical values of the masses and dominant 
branching ratios of the low mass spectrum of the heo 



Range 

Mass [GeV] 

BR1 (%) 

BR2 (%) 

ai 

[9.9,10.4] 

10.0 

T + T~ (81) 

99 (16) 

hi 

[53,59] 

55.7 

aidi (72) 

bb (23) 

ll2 

[122,123] 

122.8 

bb (65) 

WW* (17) 

Xi 

[54,60] 

57.8 

- 

- 

X2 

[57,76] 

69.8 

Xi«i(75) 

Xi (25) 

X3 

[107,136] 

121.9 

Xihi (80) 

X 2 z* (10) 

X4 

[166,209] 

179.5 

X2^ (80) 

X1^ (14) 

X5 

[225,254] 

236.6 

xtw (60) 

vv,U (38) 

xi 

[104,132] 

116.3 

X 2 IV* (78) 

XiW* (22) 

xi 

[225,255] 

237.1 

tv, tv (40) 

X 3 W (38) 

ti 

[313,391] 

335.6 

Xtb (75) 

X3t (15) 


Table 3. Masses and the two dominant decays of the lighter 
part of the heo benchmark point spectrum obtained by 
NMSSMTools4. The column titled “Range” is the mass range 
of the points surviving all scan constraints. 


benchmark, see Table [3] The ai and hi of the bench¬ 
mark avoids the LHC search exclusion for straightforward 
reasons. Both ATLAS and CMS have searched for gluon 
fusion gg —» a —> but critically omit the T re¬ 

gion and therefore cannot exclude m ai ~ 10 GeV (5SJ 
[60] . ATLAS has searched for gluon fusion gg —► h —> aa 
for 2 m T < m a < 2 ms but does not report limits for 
m,h < 100 GeV [JO]. CMS has searched for the same 
channel, but only reports limits for nih > 90 GeV with 
m a < 2m T IMi or for the hi 25 with 4 < m a < 8 GeV 
m- More decisively, the gluon fusion cross sections for ai 
and hi production in the benchmark are greatly reduced 
relative to the hi25- 

In the neutralino and chargino sector, both ATLAS 
and CMS have studied X 2 X 1 " production [5511551155] . For 
example, the searches which assume decays to sleptons or 
to bosons also assume that m x + = m X2 , motivated by 

models with a bino-like xi and wino-like X 2 and Xi~- But 
in heo the xi is singlino, and manifestly m x + ^ m X2 . The 

X 2 X 1 " searches which assume dominant decays to sleptons 
cannot exclude he 0 where > m x +,m X2 . Such searches 

might be sensitive to X 5 X 2 events, but here the cross sec¬ 
tion is reduced and the final states are more complex. Of 
the searches which assume dominant decays to bosons, 
only the W\iZ\i final state case applies. In this case 
both W and Z are very far off mass shell in heo , in which 
case it is unlikely that the exclusion can apply. 

In the stop sector, both ATLAS and CMS report ex¬ 
clusion. For a summary of the ATLAS results, see [571155] , 
For a bibliography of CMS results see [69] , No exclusion 
is given for the NMSSM, however, and exclusion for sim¬ 
plified models cannot be easily interpreted in the NMSSM 
context. For example, the analyses which assume t —> t\ 1 
with 100% branching ratio cannot exclude heo, fc> r which 
this branching ratio is 0(10 3 )- heo does contain t -> fx3 
with a branching ratio 0(1O -1 ), but the subsequent X 3 
decay produces a much more complex final state with less 
missing energy than assumed by the searches. The stop 
pair searches which assume t —> b'xf with 100% branching 
ratio assume very specific cases of mass relationships be- 
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Sample 

T max 

Analysis 

1 dtC [fb L ] 

XX 

0.6 

atlas_conf_2013_035 

20.7 

t\ti 

0.5 

atlas_conf_2013_061 

20.1 

99 

12.5 

atlas conf 2013 061 

20.1 

XX 

0.1 

cms_1303_2985 

11.7 

titi 

0.6 

cms_1502_06031 

19.5 

99 

1.2 

cms_1303_2985 

11.7 


Table 4. Maximum exclusion rmax determined by Checkmatel 
of the analyses most sensitive to /i 60 ) of all the validated AT¬ 
LAS and CMS Run 1 analyses. Note that neither stop nor 
chargino/neutralino pair production is ruled out for heo- 


tween stops, neutralino and charginos which do not hold in 
heo- Moreover they assume xt _■ > Wxi- But the dominant 
branching in heo is Xi kF *%2 — > W*Z*x i, producing 
less missing energy, and two gauge bosons which are very 
far off mass shell in comparison to the search assumptions. 

In order to evaluate quantitatively the heo exclusion at 
LHC Run 1, wc run all 16 (3) presently validated ATLAS 
(CMS) analyses in Checkmate 1.2.2 [701171117211 73 II 74 ll75i on 
generated he q events. Event simulation of the gluino, stop 
and chargino/neutralino pair production is carried out 


From Table m the decays h -a xt b -. xt X 3 W, 
X3 "A Xi^i) and hi ► 2oi proceed with branching ratios 
of 75%, 38%, 80%, and 72% respectively, while t\ -A X3 b 
proceeds with branching ratio of 15%. This makes stop 
pair production with t-\ -A xt b -A X3^ b or t-\ -A X'it and 
X3 Xi^i * Xi2ai promising channels for discovery if 
the a\ can be successfully reconstructed, for example in 
the relatively rare but very clean a\ —> g + channel. In 
heo this decay proceeds with a branching ratio of 0.3%. 
These remarks also apply to stop pairs produced in gluino 
pair production with g —> tit. 

In heo stop pair production, the cascade dominantly 
contains two top quarks. In gluino pair production it dom¬ 
inantly contains four top quarks. These are strong handles 
on any potential background. Some top pair ti background 
may be irreducible, but other backgrounds should be neg¬ 
ligible. 

We now describe a targeted study of the sensitivity to 
heo at the LHC. Signal events are generated with Pytliia 8 
as described in Section [5] Background ti events are also 
generated in Pythia 8 . Fast detector simulation is performed 
with Delphe3.2.0 [71] . The Delphes3 detector card for CMS 
is modified to reproduce the tight electron, tight muon 
and b tag efficiencies reported by CMS ,79 ,.80/81]. The 


with Pvt hi a,8.205 [761(771 . The SLHA file produced bv NMSSMignal selection seeks the decay on -A g+p in gluino 


Tools4 for heo is used with Pythia 8 , which features a ded¬ 
icated NMSSM model with functionality for SLHA input. 

See Table 0] for the exclusion r max , the ratio of the 
95% confidence level lower limit on the he o signal presence 
to the measured 95% confidence level limit, of the analy¬ 
ses with maximum sensitivity to heo chargino/neutralino, 
stop and gluino pair production. Only for gluino pair pro¬ 
duction is r max > 1, indicating that both ATLAS and 
CMS have ruled out a gluino with nig a: 611 GeV in heo 
but neither has ruled out the stop and chargino/neutralino 
sectors of heo- However, since nig is determined by the 
gaugino mass M 3 , which can be easily increased without 
otherwise impacting the lower energy he 0 phenomenology, 
we simply assume nig « 855 GeV or greater since this re¬ 
duces the gluino pair production cross section by a factor 
of 13 relative to heo- 

Note that r ma x = 0.6 for cms_1502_06031, which ex¬ 
hibits a 3(7 excess in the low dilepton mass region [78]. If 
the heo stop mass is reduced such that the stop pair pro¬ 
duction cross section is enhanced by a factor of 1.5, then 
this analysis becomes sensitive to heo with the reduced 
to,-« 315 GeV. 


6 Collider Signature 

Since the stop is relatively light in heo , the cross section 
for pair production is large and makes cascade produc¬ 
tion of the a\ and hi accessible at the LHC. Gluon fusion 
production of ai and hi is less promising. The reduced 
tthi (tta 1 ) coupling, which appears in the gluon fusion 
top loop, is of order 0 ( 1 O -1 ) ( 0 ( 1 O -5 )) relative to the 
SM ttHsM coupling for a SM Higgs boson of the same 
mass. 


and stop pair events and employs a standard selection 
for semileptonic top pair events, together with a selec¬ 
tion for ai —» p + p~, in which one top quark decays via 
t -A bW -A Mv and the other via t -A bW -A bqq'. The 
requirements for the Run 1 analysis are these: 


— exactly one tight electron with Et > 25 GeV and no 
isolation requirement 

— missing transverse energy E™ lss > 85 GeV 

— four or more jets with Et > 20 GeV, at least two of 
which are 6-tagged 

— two or more tight muons with pt > 2 GeV, no isolation 
requirement and do/<Jd 0 < 5 

— zero net charge and 9.7 < m^+^~ < 10.3 GeV in the 
leading and subleading muons 

The ai candidate is then reconstructed from the leading 
and subleading muons. The muon azimuthal impact pa¬ 
rameter significance requirement do/<Jd 0 < 5 ensures that 
the muons are consistent with prompt production. 

For the Run 2 analysis, we assume y/s = 14 TeV and 
f dtC = 300 fb -1 . We use the Delplies3 simulation with 
mean pileup 50. The selection is identical to the Run 1 
analysis except that the electron, jet, and muon thresholds 
are raised to 30 GeV, 30 GeV and 4 GeV, respectively. 

After full signal selection, the SM top background is 
nearly negligible. Multiple jet events produced by QCD 
have not been simulated, but with the nominal selection 
this background is expected to be very small. In data, the 
nonpeaking heo events in the candidate ai distribution can 
be mistaken for QCD multijets events, however, so these 
are considered background in the significance calculation. 
In the SM top background, the candidate ai muons orig¬ 
inate from r lepton, D meson or B meson decays. In the 
nonpeaking heo events, they originate either from SM r, 
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Process 

Run 

1: C = 5 fb 

=T 

Run 1: 

£ = 20 fb" 1 

Run 2: 

O 

O 

CO 

II 

fb" 1 


<T 7T [pb] 

N n 

N p 

08T [pb] 

N n 

N p 

(T14T [pb] 

N n 

N p 

SM tt 

173.6 

1.6 

0 

247.7 

8.5 

0 

966.0 

286.0 

0 

heo titi 

0.7 

2 .6(0.8) 

2.2(2.5) 

1.1 

20.1 

19.2 

5.9 

892.5 

647.9 

S/y/B 

1 .0(1.6) 

3.6 

18.9 


Table 5. NLO cross sections for t\t\ and ti production at the LHC in Runs 1 (y/s = 7,8 TeV) and 2 (y/s = 14 TeV). We 
omit the gg yields. Also shown are the expected yields for peaking events (N p ), yields for nonpeaking events (N„) and signal 
significances after the full signal selection described in the text. For y/s = 7 TeV, we show in parentheses yields and significance 
for a variation of heo in which the decays outside of the effective tracking volume. 


B or D decay or from NMSSM a\ —> t m t, y 2 —> Xi 
X + -A XU", or X + gu. 

The proportion of peaking to nonpeaking signal events 
is sensitive to the details of heo- For example, if the slep- 
ton masses are raised above the threshold for decay from 
X 2 : then the peaking signal is enhanced and the nonpeak¬ 
ing signal is reduced. Similarly, if the branching ratio for 
X 2 —> Xi a i is raised at the expense of X 2 —> XiZ*, the 
peaking signal is enhanced and the nonpeaking signal is 
reduced. Finally, if the xi width is sufficiently small, its 
decay vertices may lie outside the effective tracking vol¬ 
ume, making nonpeaking background from X 2 effectively 
invisible. 

Pythia8 is a leading order generator, but next to lead¬ 
ing order cross sections obtained by the LHC SUSY Work¬ 
ing Group |B2II53| are used to normalize the event yields. 
See Figure[5]for the reconstructed a± mass distribution af¬ 
ter full signal selection, where the distribution for a varia¬ 
tion of heo in which the X 2 decays outside of the tracking 
volume is also shown. See Table [5] for expected peaking 
and nonpeaking event yields and signal significances after 



Candidate a, Mass [GeV] 


S 

Fig. 2. Reconstructed candidate a 1 — > g + g mass distribu¬ 
tion after full signal selection assuming heo (x3 for titi and 
/13 for gg) at the LHC for y/s = 7 TeV and J dttC = 4.7 fb" 1 . 
Also shown (dashed) is a variation of heo in which the X2 de¬ 
cays outside of the effective tracking volume. The fits employ 
a Gaussian signal model and a polynomial background model. 


full selection for Runs 1 and 2 at the LHC. A targeted heo 
signal selection yields sensitivity even at the LHC Run 1. 

The advantages of the International Linear Collider 
(ILC) for studying low mass NMSSM Higgs bosons has 
been noted in m- At the ILC, the standard Higgstrahlung 
production channel e + e" —> Zh\ is suppressed in the 
NMSSM due to the measured SM-like /1125 —> ZZ* signal 
strength and the NMSSM coupling sum rule XTi=i £zzh = 
1 [14] . Instead, in heo we note the possibility of resonant 
production e + e" —> a\h\, with cross section of several 
hundred picobarns at y/s = rriz- For y/s = 500 GeV, pair 
production of all neutralinos and all charginos is accessible 
with cross sections nearing a picobarn, as well as a\h\ and 
Zh 2 production cross sections of about a hundred femto- 
barns. 

In heo the Za\h\ coupling is small enough to have 
evaded LEP searches [MlfRol but large enough to be pro¬ 
duced copiously at the ILC running on the Z pole. The 
ILC sensitivity to heo in operating scenarios described in 
m defined by beam polarization, luminosity and y/s and 
will be evaluated in a forthcoming companion study. 


7 Conclusion 

We have reviewed the motivation for a natural NMSSM 
with a slightly broken PQ symmetry and a lightly coupled 
singlet featuring a light singlet pseudoscalar ai, light sin¬ 
glet scalar hi and a light singlino LSP xi DM candidate 
annihilating via X 1 X 1 =► bb. 

A random parameter space scan is performed subject 
to a full suite of experimental constraints, including the 
anomalous muon magnetic moment, the DM relic density 
and collider searches. Surviving points are characterized 
by low fine tuning, and abundant pseudoscalar a 1 pro¬ 
duction identifies the benchmark point heo■ In addition 
this benchmark features light stops and light higgsinos, 
all characteristic of Natural SUSY. 

The benchmark avoids the current LHC exclusion lim¬ 
its. For the a\ and hi, this is due to the reduced gluon 
fusion cross sections. For other SUSY searches, this is pri¬ 
marily due to the search assumption that stops, neutrali¬ 
nos and charginos will decay directly to the LSP xi with 
no intermediate SUSY particles in the decay chain. But in 
the heo benchmark the xi is singlino and couples weakly 
to the rest of SUSY. Thus due to the light mass spectrum 
the decay chains can contain many intermediate SUSY 
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particles, making the final states more complex with less 
missing energy than in the simplified search scenarios. 

Finally, we report that the potentially fruitful discov¬ 
ery channels at the LHC for the benchmark considered are 
stop and gluino pair production with either t\ —> xt & ~i► 
X3Wb or t\ —> an d X3 -A Xift-i — > xi a i a i- We con¬ 

clude with a fast simulation study that with a targeted 
signal selection the LHC may already be sensitive to heo 
in Run 1. We have also pointed out the possibility to ob¬ 
serve at the ILC resonant e + e _ —> a\h\ at yfs = mz 
and pair production of all neutralinos and charginos at 
y/s = 500 GeV. 
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